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for the transitions 7*7rp and 7*7r°7. We discuss the role of the sub-dominant co- 
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CS ' 1. QCD Modeling of Mesons 

Q^ i The Dyson-Schwinger equation (DSE) approach [1|] to non-perturbative QCD 

0^ I modeling of hadrons and their interactions ||] combines truncated QCD equations 

r^ I of motion for propagators and vertices with infrared phenomenology fitted to a few 

^ i| key low-energy observables. Parameter-free predictions for other hadron properties 
and observables are then used to develop and test our understanding of hadron 

Q I physics at the quark-gluon level. The DSE for the fully-dressed and renormalized 

<^ . quark propagator in Euclidean metric is 



S-\p)^Z2[ij-p + mo{A)]+Zi^J Al.g^D^,{j>~k)-i,S{k)Vl{k,p), (1) 

H ; 

5^ , where too(A) is the bare mass parameter and A characterizes the regularization 

mass scale. In the DSE approach the renormalized dressed gluon propagator Df^^{q) 
and dressed quark-gluon vertex T^{k,p) are constrained in the UV by perturbative 
results and are representated by phenomenological IR forms with parameters fitted 
to selected pion and kaon observables. 

Mesons are generated as bound states of a quark of flavor /i and an antiquark 
of flavor /2 via the Bethe-Salpeter (BS) equation 

r(p; P)^ J ^K{p, q; P)Sf, {q + ^P)r{q; P)Sf., {q - ^P) , (2) 

where ^ -f f = 1 describes momentum sharing. The kernel K is the renormalized, 
amputated qq scattering kernel that is irreducible with respect to a pair of qq lines. 



The present stage of QCD modeling truncates K{p, q; P) at the ladder approxi- 
mation and couples quark color currents with bare vertices and an effective gluon 
2-point function in Landau gauge. The latter involves acsiq^) for interpolation 
between the 1-loop pQCD result in the UV and a phenonienological enhancement 
in the IR. The treatment of the quark DSE that is dynamically matched to this is 
the bare vertex or rainbow approximation; the axial vector Ward-Takahashi iden- 
tity is then preserved, Goldstone's theorem is manifest, and realistic pion and kaon 
solutions are ensured |^. 

The general form of the pion BS amplitude is 

ti{k-P)=T^^^[iE^[k-P)+rF^[k-P)+^k-PG^[k-P)+a^,k^P,H^{k-P)], (3) 

and the first three terms are significant in realistic solutions 0. With S^^{p) = 
i}tA{k^) -I- B{k^) -I- TO, the axial Ward-Takahashi identity in the chiral limit specifies 
the dominant amplitude as Ejr{k;P = 0) = BQ{k'^)/f^; the ampfitudes F^ and G,^ 
are related to A{k'^) and its derivative in a less direct manner H|. 

In Euclidean metric, the BS meson mass-shell condition requires that S{p) be 
evaluated in a certain domain of complex p^ . The required domain for meson decays 
and form factors can be quite demanding. To facilitate such studies, we make use of 
an analytic parametrization of the numerical solutions of the quark DSE to represent 
S{p) as an entire function in the complex p^-plane describing absolutely confined [g| 
dressed quarks. Typically five parameters are used to achieve a good description 
of pion and kaon observables: fyr/K'^ '^Tr/if! (99); the tt-tt scattering lengths; the 
charge radii r-jr, rj^o and r^±; and the pion charge form factor |^,0. Current efforts 
along this line are concentrated upon vector and axial vector mesons and hadronic 
and semi-leptonic decays. Here we outline recent studies of several electromagnetic 
transitions and also the pirir coupling. 

2. The pirn Coupling Constant 
The first term in a skeleton graph expansion of the ptttt vertex y,^ is 

AAP,Q) = -P.9prrnF{Q^) (4) 

= 2iVetr,y'^r^(fc"';-P+)5(g')re(fc';g)5(<z")r.(fc";F_)5(g"'), 

where P± ^ P ± Q/2 and both pions arc on the mass-shell. The first argument of 
each BS amplitude is a relative qq momentum (we choose equal partitioning) and 
the second argument is the incoming meson momentum. After the loop momentum 
k is specified in terms of one internal momentum, the others are easily deduced. By 
definition we have F{Q^ = —m?) — 1. 

For this ptttt study, we employ approximate F^ and T^, obtained from a rank- 
2 separable ansatz H for the ladder/rainbow kernel of the DSE and BSE. The it 
is properly massless in the chiral limit as required by Goldstone's theorem; but 



the full consequences of the axial Ward-Takahashi identity M are not preserved. 
Parameters are fit to mjr/K and Jt^/k- The results have the form |^ 

r,(fc; P) = Z75 /(fc^) A^ - 75 f I{k^) AJ , (5) 

re(fc; P) = fcj g(fc2)Aj + i7j /(fc2)A^ + ^75e^,;,^ ^^k^P^ g{k^)\^ . (6) 

Here g{k'^) = [A{k'^) — V\/a and /(A:-^) = B{k'^)/b with a and 6 given in Ref. Q from 
the quark DSE solution. The relative strength of the Ai is given by the separable 
BSE solution which requires a phenomenological UV suppression. In contrast to 
previous work H, here we normalize the A^ in the canonical way [H without such 
a suppression in any of the momentum dependent quantities. The results for gp-m^ 



Table 1: 



Qp-K-K 



calculation and contributions from meson covariants. 



gp^^ = 7.32 [expt 6.05] | 


TT Covariants 


p Covariants 


75 184% 
757 • Q -84% 


7m 98.5% 
75eM ikQ 1.4% 
/Cm 0.1% 



are given in Table |l| where an error in our previous work |10| is corrected. The 
empirical value associated with the p —>■ tttt decay width is overestimated by 21%. 
The integral for gpT^r is significantly influenced by the normalization of r^r ; the form 
of the separable amplitude Eq. (g|) may be too primitive to generate a realistic norm. 
Improved BS amplitudes |^,Q are available for future work. The pseudovector n 
component is much more important here (-84%) than for it is for ttt-tt and /,r (25- 
35%). The sub-dominant p amplitudes make only a minor correction. 

3. The 7*7r/9 Form Factor 

The isoscalar 7*7rp meson-exchange current contributes significantly to electron 
scattering from light nuclei. Our understanding of the deuteron EM structure func- 
tions for Q^ ~ 2 — 6 GeV^ requires knowledge of this form factor |1J] . The general 
form of the vertex, and the explicit quark loop that arises in the impulse approxi- 
mation generalized to include the dressing of the propagators and vertices, is 



Ap,.(P, Q) = -i ep^a/3 Pa( 



2Nc 



trc 



(2^ 



f,(fc" 



!57-P^(Q') (7) 

-P+) 5(g') r.(fc'; Q) S{q") r^(fc"; P_) S{q"') . 



The momentum notation is the same as for the previous pmr case. The dressed 
photon-quark vertex is taken to be the Ball-Chiu [|l^ ansatz 



r.(fc;Q) = -^(^+ + A_)+^ 



ifi{A^ -A+) + B^ -B^ 



(8) 



where f± = f{k±) and k± — k ± ^. This form obeys the Ward-Takahashi identity 
(WTI) and the relevant symmetries and is conveniently determined completely in 
terms of the quark propagator. It then follows that Qi, A^i, = 0; the ^np current 
is conserved. With the above separable model tt and p BS amplitudes, along with 
the associated quark propagator parameterization |^, we obtain g^T^p = 0.45 in 
reasonable agreement with the empirical value g°^p — 0.54 ± 0.03 from p decay. 
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Figure 1: The 'y*7Tp transition form factor with vr and p on shell. 

The pseudovector tt amplitudes i<V and Gt^ defined in Eq. (0) generate the 
correct asymptotic behavior of the pion charge form factor |1J]. Also their UV 
relationship k'^G^^ik'^) — > 2Fjr{k'^) implements convergence for the /^r integral (l4|. It 
is impossible to generate these properties in the previous separable ansatz approach; 
we use here the approximating forms |14| of a numerical solution Q 



E^k'') 






F^k^ 



E.{k^) 

iiOA 



G.(fc' 



2F^{P) 

fc2 + M^y 



iJ,(fc2)«0, (9) 



where N^^ is fixed by the standard BS normalization condition M. We use the 
quark propagator parameterization (set A) |14| associated with Eq. (S). There is 
no dynamically matched F^ available so we simply take Eq. (0). We then obtain 
g-^Trp — 0.708. Since the canonical p covariant is so dominant, we expect that the 
mismatch is simply one of normalization and that the produced F{Q'^) should be 
quite realistic. The result in Fig. n^is much softer than the vector meson dominance 
(VDM) prediction. The available data for elastic EM deuteron form factors A{Q^) 



and B{Q'^) in the range 2 
previous 7*7rp vertex result | 



6 GeV has been shown [IJ] to strongly favor our 
] . The present work is less phenomenological, employs 



more realistic representations of the tt and p, and produces a harder form factor. 
The influence on the deuteron form factors remains to be determined. 
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Figure 2: The 7*7r°7 transition form factor. The data are taken from Ref. ^ 
(CELLO) and Ref. jl§ (CLEO). 



4- The 7*7?*^ -^ 7 Transition 

The coupling constant for the tt*^ -^ 77 decay is given by the axial anomaly and 
its value is a consequence of only gauge invariance and chiral symmetry in quantum 
field theory. The form factor of this anomalous transition is not dictated by symme- 
tries; it is of interest as a test of our ability to model nonperturbative QCD because 
of the relatively simple hadronic dynamics that is involved. In the asymptotic UV 
region, one expects a simple result dictated by the known electromagnetic coupling 
to current quarks and the intrinsic properties of the pion. 

The general form of the vertex allowed by CPT symmetry, and the explicit quark 
loop that arises in the impulse approximation generalized to include the dressing of 
the propagators and vertices, are 



^liuiP,Q) ^ i—^ <^tiual3 PaQl3 gTTjf F{Q'^) 



(10) 



^tr. 
3 



(2^ 



S{q')r,{k'; Q) S{q")T,{k"; -P - Q) S{q"')TAk"' ; P). 



The 7* momentum is Q, and the other photon and the pion are on the mass- 
shell. With the convenient choice k'" = k, the other internal momenta can be 
determined. The Ball-Chiu ansatz Eq. (pi) is used for the dressed-quark-photon 



vertices Ti, and F^. The chiral limit anomalous decay 7r° — > 77 gives g^~.^ = 1/2 
providing an excellent account of the 7.7 eV width. The form factor defined by 
Eq. (nG) satisfies F{0) = 1; it should not be confused with a different quantity 
F{Q^) = F{Q'^)/'^T^'^ fir which contains the non-tensor strength of the transition 



matrix element M„ 



2Kr.i, and in terms of which the CLEO E6(l data and some 



theoretical works are expressed. 
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Figure 3: The 7*7r*'7 transition form factor times Q^. 
are described in the text. 



The asymptotic limits marked 



Here we update an earlier study |l5|] by using the more realistic Fjr in Eq. (||) 
along with the associated quark propagator (set A) (l4|. We obtain .97177 = 0.4996 
at the physical 771,^ value. The form factor is displayed in Fig. |[ Although the 
anomalous coupling strength is correct, the calculation does not fall off fast enough 
in the infrared and significantly overestimates the data. It is unlikely that our 
description of the pion is responsible for this; the sub-dominant amplitudes i^^r, Gtt 
give a negligible contribution here while the closely related dynamical quantities, r,r 
and /tt, are well described [Q . The removal of dressing at the virtual photon vertex 
produces the dashed line. The better agreement with the data is rather fortuitous 
since the coupling constant is reduced to 70% of the former correct value. The pion 
transition radius from the impulse approximation (~ 0.48 fm) is clearly less than 
that suggested by the data ('^ 0.65 fm); a similar underestimate also occurs for the 
pion charge radius in this approach. 

Our results for Q^F{Q^) are displayed in Fig. 0. The approach to the asymp- 
totic limit of the present dressed quark loop Eq. dTO) is quite slow and is governed 
by the following considerations. After the mass-shell conditions for the pion and 
the real photon are realized, one finds that in the domain fc^ < 1 GeV^ ^ Q^ of 



integral support dictated by Tj^ik'^), the leading behavior of the momenta for each 
quark propagator is: (g')^ {q"T ^ 0{kQ), and {q")^ = ^[l + 0{k/Q)]. Using 
A{p'^) = 1 + 0{l/p'^) and B{p'^) - ©(l/p^), one finds from Eq. (||) that the lead- 
ing behavior of each of the photon vertices is T^ = ~ij^[l + 0{l/kQ)]. The pho- 
ton insertions of the diagram effectively collapse to a point axial vector with the 
ri,S{q")rf^ leg of the loop having the leading form 

^ 7^ ^ 7m = Q2 '^M'-a/J IblaQp H , (11) 

where the current mass is ignored and the terms not shown do not survive the spin 
trace. The loop integral coupling the tt to 757a gives exactly f-„Pa/2> by defini- 
tion [Q|. This gives the asymptotic hmit k^i,{Q'^) — > ^ i^^ujap PaQp[^ + 0{1/Q)] 
or the form factor limit ^t:^ fl + 0{l/Q), in agreement with Refs. |l|,|§. The 
limit marked by Stt^/^ is from pQCD factorization fl^. It is evident that the sub- 
dominant amplitudes Ft^ and G-^ make a minor contribution although the integrated 
effect via the produced /^ value is some 30%. 

With the hard leg Ti,S{q")Vfj_ taken to be bare the result is the dot-dashed 
line in Fig. y. This is consistent with the evident sub-leading correction in the 
propagator denominator being (k/Q) ^ 30% at Q^ ~ 10 GeV^. The dressing of the 
hard propagator S{q") contributes little to the dot-dashed curve; thus the difference 
between the three relevant curves illustrates the persistent contribution from photon 
vertex dressing. Since the Ball-Chiu vertex is exact at both Q^ — and the UV 
limit, and only the longitudinal component is correct for all Q^, it is possibly the 
deficiencies of this ansatz at infrared and intermediate momenta that are being 
exposed in the present study. This is also the preliminary finding from a study of 
the ladder Bethe-Salpeter solution for the vector vertex pl| . 

5. Summary 

The results presented here suggest that the present approach to modeling low- 
energy QCD can capture the mechanisms that dominate infrared physics. The 
sub-dominant pseudovector terms in the pion BS amplitude make a much stronger 
contribution to gpTr-K than they do to ttitt, Jtt', the contributions to the 7*7r°7 and 
7*7rp transition form factors are minor at low and intermediate momenta. The new 
result here for the 7*7rp form factor should be applied to the deuteron electromag- 
netic form factors where this meson exchange current is still a serious ambiguity. 
Our examination of the pion axial anomaly form factor with the Ball-Chiu ansatz 
for the dressed photon-quark vertex clarifies the slow approach to asymptotic be- 
havior and suggests that deficiencies in this ansatz at low momenta may be evident. 
A study of electromagnetic radii should be made from this perspective. 
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